Tensor of coherences parameterization of multiqubit density operators for 

entanglement characterization 



00 

o 



Oh! 



PACS numbers: 03.67.Mn, 03.67.-a, 03.65.Ud 



Claudio Altafini 
SISSA-ISAS 
International School for Advanced Studies 
via Beirut 2-4, 34014 Trieste, Italy Q 

For multiqubit densities, the tensor of coherences (or Stokes tensor) is a real parameterization 
obtained by the juxtaposition of the affine Bloch vectors of each qubit. While it maintains the 
tensorial structure of the underlying space, it highlights the pattern of correlations, both classical 
and quantum, between the subsystems and, due to the affine parameterization, it contains in its 
components all reduced densities of all orders. The main purpose of our use of this formalism is to 
deal with entanglement. For example, the detection of bipartite entanglement is straightforward, 
' as it is the synthesis of densities having positive partial transposes between desired qubits. In 

' addition, finding explicit mixtures for families of separable states becomes a feasible issue for few 

, qubit symmetric densities (we compute it for Werner states) and, more important, it provides some 

insight on the possible origin of entanglement for such densities. 

OJQ. 
< 

^ '. I. INTRODUCTION 

o\ 

The main purpose of this work is to discuss a parameterization of a multiqubit (pure or mixed) density operator and 
' to show its usefulness in dealing with classical and quantum correlations. The principle behind the parameterization 
is the same of the so-called vector of coherences of widespread use in modeling TV-level density operators, Q, H, |H If^ 
2^. It consists in choosing a complete orthogonal set of Hermitian operators and in considering the corresponding 
. real vector of expectation values in place of the density matrix. Here we use the same idea, but respecting the 
tensorial structure of the density, hence working with a (real) tensor of coherences. Each qubit is parameterized as 
an affine Bloch vector, and the tensor is just the juxtaposition of affine Bloch vectors. Maintaining the tensorial 
structure has several advantages, for example it makes the pattern of the "total" correlation between subsystems 
' totally straightforward to see. For "total" correlation we mean both the classical and the quantum ones. As a matter 
^ ■ of fact, we will see that the correlation between subsystems is encoded in the terms of the tensor. Crucial to the 
^ understanding of this point is the role of the affine component and of how it enters into the compounding of the 
different qubits. In fact, as we use homogeneous coordinates to deal with the affine term, the key simplification is that 
^ ' tracing over one of the qubits simply corresponds, up to a scale factor, to choosing the "0" (i.e., affine) component for 
• i-H , the corresponding index. Hence, because of the affine parameterization, reduced densities are naturally represented by 
' means of the tensor of coherences parameterization and the tensor itself consists of the entire hierarchy of correlations. 
H I Also the scale factor has a natural interpretation: it corrects the trace norm of the completely random state when 
■ - - I passing from a density to a reduced density. Further advantages are in the simplicity of the geometric picture for 
multipartite systems and in the possibility of using multilinear algebra ideas in a more strairfitforward manner. 

Obviously the parameterization we consider is not really new; it was treated in detail in [lOj, used extensively for 
exainple by Mahler and co-workers [23, 13 (where the basis elements we use are referred to as cluster operators) or in 
[Tsl ll^ (where it is referred to as the Stokes tensor) and more or less implicitly in many other papers, cf. hj for 
related material. For example in the NMR literature Q it goes under the name of product of operators basis. What 
is new is its use in understanding multiparty entanglement, see [M ll2l[T9.. 26j for an overview of research in this field. 

The simplest (to detect) type of entanglement is bipartite entanglement, for which there exist a necessary and 
sufficient condition, the so-called positive partial transpose (PPT) criterion of 0,|23|. In the tensor of coherences 
parameterization, the PPT criterion has a very simple formulation and, more important, it becomes completely trivial 
to construct densities satisfying PPT between all pairs of subsystems. Hence one can focus on the class of entangled 
PPT densities, which are characterized by the more subtle bound entanglement |3- On the other hand, also the 
construction of systems having certain patterns of bipartite entanglement not satisfying the PPT criterion (NPT 
entanglement) is rather simple. 
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If the total correlation, classical plus quantum, is directly depicted in the tensor of coherences, the distinction 
between the two types of correlation remains however an elusive issue, although as we will see in the examples, the 
parameterization allows to suggest what is happening in an entangled state. Consider a one-parameter family of 
densities ending in the maximally mixed state. Close to such extreme the state is certainly separable ^2^. Using 
the tensor of coherences it is not too difficult to construct an explicit (one parameter) convex mixture of product 
states for it. In all the examples of entangled families we have tested, the convex combination found is such that 
it induces cancellations between the corresponding mixtures of reduced densities. If close enough to the complete 
mixing these cancellations are harmful, far from it it may happen that the one parameter density is well-defined while 
some of the reduced densities (which, again, being canceled do not explicitly appear unless one wants to construct 
the mixture explicitly) are not anymore well-defined densities in the sense that their trace norm is too big (some of 
its eigenvalues become greater than 1). Hence one source of entanglement is that not well-defined components give 
rise to a well-defined compound system. The complication is obviously that due to the nonuniqueness of the mixture 
representing a given density, it is an hard problem to exclude that any other convex combination will suffer from the 
same "unfeasibility" problem. For low rank systems, the tensorial notation helps in finding such convex combinations. 
For example we could easily compute a mixture for the Werner states valid in the whole separability interval. 

The parameterization into vector (and tensor) of coherences is natural only for qubits. For a fc-level system, in 
fact^he vector of coherences of the density operator is not free to evolve on the corresponding (affine) ball in R*^ 
sec 5, 18] for hints on this point. Of course qubits are by far the most popular systems in quantum information 
processing. 

One may argue that the dimension of the state tensor grows as 2^" with the number n of qubits and hence that 
expanding densities explicitly into a complete basis becomes rapidly cumbersome. The exponential growth of the 
number of degrees of freedom available concerns however all densities, regardless of the representation used. While 
this fact is immediately evident using our notations, it may go unnoticed using some standard parameterization. Of 
course in a problem like detecting entanglement all degrees of freedom of the state may come into play, therefore we 
find it convenient to have them all explicitly expressed. 

In next Section the tensor of coherences is introduced and correlations, both classical and quantum, are discussed in 
its terms. For sake of notation simplicity, we consider in some detail the geometry of the 2-qubit case. The extension 
to n-qubit densities is straightforward. In Section lllll several examples are treated. We construct explicit mixtures 
for Werner states and for a tripartite family ending into the bound entangled state of 3. An example of how to 
construct (and analyze) NPT tripartite entanglement is also proposed. 



II. TENSOR OF COHERENCES PARAMETERIZATION FOR p 

Given n qubits living on the Hilbert space (7i^)®" of dimension 2", the corresponding density operator is a 2" x 2" 
positive semidefinite Hermitian matrix p such that tr (p) = 1 and it has 2^" — 1 degrees of freedom. We construct 
for p a basis borrowed from the literature on NMR spectroscop y w here it is normally referred to as the product of 
operator basis 0. Similar bases are discussed for example in [Tolll5lll7ll2^ . In terms of this basis, studying densities 
is equivalent to studying tensors of directly observable real parameters. 

A word on the notation: we use the symbol "p" for density matrices and "g-'" (possibly with a multiindex) for the 
components of the tensor of coherences. The superindex is always a tensor (multi)index; for powers of we use an 
extra round bracket. For the tensor, we also use the summation convention over repeated indexes, always in the range 
{0,...,3}. 



A. One qubit 



The rescaled Pauli matrices A, 



Ai 



^aj, j = 1, 2, 3, 
V2 [l 0. 



1 

V2 



~i 

1 



1 

V2 



1 

-1 



plus the rescaled identity operator Aq — -^12x2 form a complete orthonormal basis (in the sense that tr (Aj Afe) = Sjk) 
for 2x2 Hermitian matrices. Fixing the trace means fixing the component along Aq. Hence p can be expressed as 
the affine 3- vector 



p = £»°Ao + g^Xi + Q^\2 + g^X^ = g>\j 
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where g' = tr (pXj), j = 1, 2, 3 and the component along Aq is = tr(pAo) = tr (p) /\/2 — 1/a/2. Since p° is 
a constant, it is normally neglected and only the Bloch vector g — {g^ g^ g^)'^ is considered. However, here it is 
convenient to keep the constant part and to represent the afhne vector in terms of a set of homogeneous coordinates, 
i.e., by means of the 4-vector g = g^ g^ Q^Y^ ■ From tr(AjAfc) = 5jk, j, /c = 0, 1, 2, 3, tr(pip2) induces an inner 
product on the parameter space (M^) given by tr(pip2) = ((gi, Q2f) = Q1Q2 + ((Pii ^2)) — ^ + {{qii 92))- The norm 
of g is then given by ||^|| — ^ {{g, g)) — ^tr (p^) — \J \ ^ II^?1P — + r^. Purity corresponds to tr (p^) — 1 i.e., 
II^pIP = i or g belonging to the sphere of radius r = call it S^y^, while the 4-vector {g^ g^ g^ g^Y belongs 

to the affine sphere (^g^ , ^i/x/l) ~ ("75' ^i/v^) ^ ^O'^Pl^t^ mixing, given by p = 5I2X2 = "^Aq, has norm 
tr (p^) = {g^Y — 5 ^i^d corresponds to ^3 = i.e., to a "sphere" of radius. All degrees of mixing are in between the 
two extremes just presented and in general the Bloch vector p e §^ for ^ r ^ . Hence we have — ^ ^ ((pi , ^2)) ^ 5 
andOs; ((ei, ^2)) s$ 1, Vpi, g2 G §^ 



B. Two qubits 



Call Ajk = Xj iS> Afc, J, k £ {0, 1, 2, 3}. Up to a normalization constant, the Ajk form the so-called product operator 
basis, see j9j, and are subdivided into 

qubit operators Aqo 

1 qubit operators Aqi, A02, A03, Aio, A20, A30 

2 qubit operators An, A12, A13, A21, A22, A23, A31, A32, A33 

Similarly to the 1-qubit case, the set of Ajk j, k e {0, 1, 2, 3} forms an orthogonal basis for all 4 x 4 Hcrmitian matrices 
(as p is now). It is still normalized i.e., such that tr (AjfeA/m) = tr(AjA; X^Xm) — tr (AjA;) tr (AfeAm) — 5ji5km for 
all j, k, I, m & {0, 1, 2, 3}. Except for Aqo, every Ajk has 2 eigenvalues each with multiplicity 2. An equivalent 
description of p is given by the 2-tensor g''^, j, A: = 0, ... 3, where g''^ = tr {pAjk) i.e., 

p = g>''Ajk = g^'^Xj ® Xk (1) 

The tensor of coherences g'^ can still be seen as the 16-vector [g^^ g^^ . . . g^^ g^^ . . . g^^Y , which is still a homogeneous 
representation of an affine 15-vector since g^^ is a constant. In fact, as Ajk is traceless for {jfc} 7^ {00} and Aqo ~ 
\^2y.2 ® l2x2 has trace 2, tr (p) = 1 implies that g^^ = tr (pAoo) = tr (p) /2 = 1/2. Again the g'^ parameterization 
lives on R^^ endowed with the Euclidean inner product one gets from the following: 



tr (p2) = tr ((^>"=A,fe)') ^ (^''')' = co^^t 1 (2) 



j,fc=0 



Following the terminology of the norm tr (p^) is a quadratic Casimir invariant of p. Following instead for example 
[23|. tr (p^) is a Tsallis entropy corresponding to the choice of parameter <; = 2 in Sq{p) = (tr (p^) — 1) /(I — g). 

Unlike the single qubit case, the subset of M}^ in which the parameters g'^ are such that g^'^Ajk is a well-defined 
density operator is not at all clear a priori and a hierarchy of nested subsets exists: 



uncorrelated C separable C entangled C "nondensity" C 



pl6 



If the density operator p is uncorrelated (i.e. p is a product state: p — pA® Pb), then also g''^ can be intended as 
the tensor product (which for scalar quantities becomes ordinary multiplication) 

g^'' = Qa ®Q% = QaA (3) 

with g-'^ describing the state of the first spin and g^ the state of the second. By "fully stretching" g''j^ ^ g%, one 
obtains still the 16-dimensional vector. Eq. is not true for correlated states: for example ^ = . . . g^^] has no 
expression of the form |j2Jl. Notice that since g\ = g'g are nonzero constants, if g^'^ = for all pairs {jk} such that 
j ^ and fc ^ 0, then p is uncorrelated. From the same argument, it follows that, even for correlated densities, it is 
always possible to write and g'^'' in the form: 

gio^gJ^^gOg^l^g^ and g""" ^ g^A ^ ■ (4) 
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The corresponding gi-^ = V^g-'^ and gg ~ v^f?*"^ are univocally determined. Therefore g^*^ = Pa ® = 5' regardless 
of the uncorrelation of p. 

The main difference with respect to other papers hke |lfl l23 | is that we include the 0-qubit and 1-qubit terms in 
the dyadic tensor structure. As a consequence, in the basis JQ), the reduced density operator is very natural to obtain 
as the partial trace operation consists simply in selecting the component of index "0" in the qubit to be traced over. 

Proposition 1 Given p = g'^hjk, the reduced density operator is given by 

PA = tvB (p) = g\\j = V2g^°Xj (5a) 
PB - tr^ (p) - e|A, - V2g'"'Xk. (5b) 

Proof. Since tr(Afe) = 6k, 

PA = trs (p) = g{Xj = gihiB (A^fe) = g^'^Xj (g) tr (Afe) = V^g'^Xj 
or — V2p-'° from |0J. Similarly for ps- □ 

In general (also for noncorrelated or nonseparated densities) the numbers tr {p\) and tr {p%) in {Tj) are partial 
quadratic Casimir invariants i.e, the quadratic Casimir invariants of the two reduced densities. The scale factor V2 on 
all components of the reduced density plays a double "normalization" role: it takes care of the trace and it modifies 
the quadratic Casimir invariant of the completely mixed state (which changes with the number of qubits). 

For an uncorrelated p, from the tensor product structure, g^'' ~ g^® g% is living on the tensor product of two 
affine spheres in M"*: 

g^^ ® g% € §L) ® K^j , < r^, rs ^ (6) 

and therefore eq. l(2Jl becomes: 



3 
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tr(p2)= [s'A®9'h) - 

\ + UAf) « + ll&ir) = + + 4) = tr {p\) tr . 

Under local orthogonal transformations, rA = const and — const. Under nonlocal rotations instead Q is not 
anymore valid. Likewise, the factorization is not valid for correlated p. 

Equation ((HJ is useful to have a geometric picture of the manifold on which g'^ lives. Expanding into its 16 
components, we have 

• g^'^ = 51^ the affine part; 

• [g g g ) ^ ^r^/^/2' 

- („m Sl2 03\ p 5^2 



• the remaining 9-dimensional vector ig^^ g^^ g^^ g^^ g'^^ g^"^ g^^ g^'^ g^^) £ §^ ^ 



With the parameterization chosen, testing correlation becomes an almost tautological issue. For example, a nec- 
essary condition for uncorrelation of p is that g^^ ^ and (p^ ^ Vp^*"' 7^ 0, j, k = 1,2, 3. In fact, assume that 
p is uncorrelated and g'^ ^ for some pair {jk}, j, fc = 1, 2, 3. If for example g'^ = then from g'^ = g^A® g% 
we must have that g^A — ^ and hence g'^ = f?^ ® f?B = 0. Another necessary condition for uncorrelation of p is that 
tr {p^a) tr {p^b ) ~ tr \p^\ That these conditions are not sufficient for uncorrelation is easily shown by examples, see 
Section ImH Following it is convenient to describe the correlation between A and B by means of a correlation 
tensor C = c^'^Ajk (homogeneous i.e, such that c^*-' — c-'" = 0) and rewrite p as 

P^PA<E)PB+C= (^g^Xj^ (g) (g'hXk) + c>^Ajk (8) 

where d'^ = (f^ — g^gg, j, k — 1, . . . ,3. Hence the necessary and sufficient condition for uncorrelation of p is that 
C = 0. On the other extreme, when the reduced densities pA and ps are completely mixed the tensor of coherences 
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contains only correlations: g^'^ = c^^ , = 1, . . . ,3. From © we can compute WqW^ = WqaW^WqbW^ + Y^j k=i {'^''^) ' 
obtaining the classical inequality 

tr < tr {p\) tr (9) 
A classically correlated or separable state is written as a linear convex combination [28j | 

s 

P^wPpA,p(E> PB.p, wP^O, '^wP = l. (10) 

p=l 

If for pure states separability is equivalent to uncorrelation, for mixed states the simultaneous presence of both classical 
and quantum correlations complicates considerably the picture. 

The convex combination reduces under partial trace operation and the statement of Proposition ^ implies 

Proposition 2 Given p = wP pA,p ® PB,p, wP ^ 0, ~ 1' reduced density operators are given by 

PA = trs (p) = vjPpA.p = V2wPg'^pXj (11a) 
PB = tTA (p) = wPpB^p = V2wPg°i^pXk. (lib) 

Proof. The result is well-known and the proof is reported here as an exercise in computing with the tensor of 
coherences. Considering p 

P = {w^Qa,iQ%a + ■■■ + w'Qa,sQ"b.s) Aoo + {w^gA,igB,i + ■■■+ w'gX^gl^,) Aqi + . . . + 
+ {w^qXi0%,i + ■■■ + WgX.gl^,) A33 

and tracing over B 

PA = ti-Bip) ^V2{{w^gXig%A + ■■■ + w"9XsQ%.s) >^o+ 

+ {w^Qa,iQ%a + ■■■ + w'gX,g%^,) \i + . . . + {w^gA,iQ%,i + ■■■ + WgX.g^J A3) 
= V2 {w' (gX.Xo + ... + gX^h) g%^, + ... + w' (gXs^o + ■■■ + qXs^s) 9%.s) 

= W^pAA + ■■■ + w'pA,s (12) 



since g% p = -i=, Vp = 1, 



□ 



It is possible to expand p — wP pA,p ® PB,p emphasizing its affine structure as follows. 

P = ^Aoo + -^wPg'^ pAjo + ■■■ + -^wPgXpAao + -^wPg^^pAoi + • • . + -^wPg% ^Aos 
+wPg\pgBpAii + ... + wPg\pg%pA33 
= ^Aoo + PA (E) -^Xo + -^Xo(S) PB +wPg\pgB.pAii + .. . + wPg\.pg%pA33. (13) 

The expression (|13|l allows to easily provide a geometric picture of the state space of p, extending the one for 
uncorrelated densities. If — wPrj^ p and — wPrs^p, then from (|13|l : 

• ^ is the affine component 

- (AO „20 „30^ — J_^, p S2 _ pap 

•[egg ) ~ ^gA(^\^/^-w h^^ ^^^ 



pg2 



rB.p/V2 



• [g g g ) - ^gB e 

. (^?" . . . ^.33) = {njPg\yB,p ■ ■ ■ wPgXpg%,p) e u^^s?,,^ ® s^,,^. 
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Notice how in the last item the "multipHcative" part of p hves on a convex sum of tensor products of spheres. Varying 
p, tb^p and , all convex combinations of tensor products of spheres of all possible radii can be obtained. Since 
any point in the Bloch ball is a good vector of coherences for a qubit and since a convex combination of compact 
convex sets is compact and convex, the "multiplicative" part of a separable p, g^^ , . . . , g"^^, lives on a compact, convex, 
infinitely generated set 29] . 

For p separable, p = pA,p ® PB,p, we have 

s 

^^{P^) = X! '^^'^''iSA.p, BA,q)){{QB,p, BB,q)) 
P,<3=1 



{wn'\\-BAA'\\QB,pf + E 



(14) 



W''w''{{QA,p, gA,q)){{gB,p, gB,q), 



{{Qa., :P : QA,q)} is an inner product between vectors of possibly different lengths: {{gA,p , QA,q}) — 5 + {{Qa,p , QA,q)) 
with gA.p G ^r^pi QA.q & ^r^,- Since ^ r^j ^ J = 9, we have by the Cauchy-Schwarz inequality 

-5 ^ -rA,prA,q ^ {{qa.p , BA,q)) ^A,pf^,9 < \ and hence ^ , £49)) < 1- 

For p separable, beside @ we also recover the necessary condition of [23 (itself a consequence of the partial disorder 
criterion, see 21]) on the traces of the reduced densities 

Proposition 3 Given p = w^pA^p ^ PB,p, ^ 0, J2p=i = 1j have that 

tr{pl) ^ tr(p2) (15a) 
tr(p|) > tr(p2). (15b) 

Proof. 

From Proposition 12 



p,q=l 



^{wn^-gA,pf+ E w^^'iBA^prBAA 



(16) 



Equation (|14() and (|16|) have the same number of terms, but each term of (|16(l is greater or equal than the corresponding 
one in lfn|) . since ||^_B,p|p ^ 1 and {{gs^p , BB.q)) ^1- □ 

Eq. p5|l could be rewritten in terms of the tensor of coherences as 

UAr^ur, UBr^ur. in) 

From in7|) . the geometric interpretation of itT^ is straightforward: UpaIP ^ imphes ^ + r\ > j+r^ i.e., the radius 
of the sphere of the compound system is bounded above by those of the reduced densities: ^ ^ min{r^, r'^} + j- 
The coefficient j takes care of the displacements of the centers of the p and pA (or pb) spheres due to the different 
affine terms. Such property is invariant to convex combinations. Notice that this is a necessary but not sufficient 
condition for separability (indeed it is not even sufficient for uncorrelatedness) of p. It becomes a necessary and 
sufhcient condition if p is pure, see [2^. 

Also the Partial Transpose operation of [23| p^^ = (T ® I2) (p) and p^^ = (I2 ® T) (p) (where T denotes the single 
qubit transposition) becomes a very straightforward operation in the chosen basis. 

Proposition 4 The partial transpositions of p become the sign changes on the corresponding index "2": 

p^' = g°%fe + ^^l'=Alfc - g2fcA2fe + ^^3%fc (18a) 
= g^%o + e^'A,i-e^'2A,2 + e^'3A,3. (18b) 

The proof is by direct computation. As is well-known |llj |. the PPT criterion provides a necessary and sufficient 
condition for separability in the two qubit case: p is separable if and only if p'^^ (and p'^^) is a density. We say in this 
case that p is PPT. Proposition ^provides the following constructive separability condition. 
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Corollary 1 A sufficient condition for PPT is that all terms of indexes '2' in the tensor are 0. 

For the tensor q^^ , the partial transposition is a linear, norm preserving operation: tr (p^) = tr ^(p^^)^^ = 

tr ^(p"^^)^^ • Hence entanglement violating PPT docs not modify the quadratic Casimir invariants of the density and 
the necessary conditions (|15|l are insensible to it. 



C. n qubits 

The notation for a density operator composed of n qubits (labeled Ai^ . . . , An) is completely analogous to the 
two-qubit case: 

p = PAi...A„ = = ^■'^■■■■'"Aji (g) . . . «) Aj„, ji, . . . , j„ = 0, 1, 2, 3 

If Ajj...j^j has k ^ n nonnuU indexes, then it is an operator on k qubits. The object g^i - -'" can still be seen as an 
n-tensor or as a 2^" — 1 affine vector having constant component g'^---^ = g^_^ ® . . . ® g^^^ — i-^=] (check that indeed 

tr(Ao...o) — tr(Ao)" = hence tr (p° - ''Ao...o) = !)■ Reduced density operators are obtained analogously to 

the bipartite case by collapsing each index being traced over to and rescaling by each time. For example, the 
reduced density operator of the k-th qubit is obtained by tracing over the other n — 1 qubits: 

PA,=tTA,...A,^,A, + ,...AAPM...Aj = = (V2r-V°-°^'^'"-%., Jfc = 0,...,3 

Similarly to the bipartite case, we have that p is said uncorrelated if 

p = PA^<S> ...(g) PA„ = Qa,^ ---^ Qa^^Ji <g • ■ • <g Aj„ = Qa.Qm ■ ■ ■ 6'A„^JiJ2...i„: 

and for p uncorrelated we have 



tr(p2) =tr{pAjti{pA,)...tT{pAj= Q + IIpaJI'^ Q + 



2 \ / 1 II :^ l|2 



QA2\\ ]---[^ + \\QA„ 



A necessary and sufficient condition for uncorrelation of p is that for all n-tuplcs of indexes {ji . . -jn}, jij ■ • ■ , jn — 
0, 1, 2, 3 

= gji^ ^ , , , ^ ^j^^ ^ ^n^ , , , ^j.^ V ji , . . . , j„ = 0, 1 , 2, 3 (19) 

Obviously the test H19|l can be restricted to n-tuples {ji ■ ■ ■ in] such that at least two indexes are nonzero. 
The density p is said separable if 

p = wPpAi,p® Pa^,p (20) 

where X]p=o '"'^ = 1, if^ ^ 0. Unlike the 2-qubit case, no neat necessary and sufficient condition for separability 
is known, although for few qubits complete classifications in various entanglement classes have been proposed, for 
example in ^ . Qualitatively, one distinguishes between entanglement due to violation of the PPT test for some of the 
indexes in {ji . . . j„} and a more subtle type of entanglement, called bound entanglement which cannot be detected 
by means of PPT. The first type of entanglement is bipartite and will be referred to as NPT entanglement. 

Proposition 5 p is NPT entangled if for some fc G {1 . . .n} the partial transpose 



^Tfc _ ^'l...jfc-lOjfc + i...j„ A . . . I ^'l---jfc-lljfc + l---jn A . 



(21) 



_ „ji---i/c-i2ifc+i...j,i A . . „. . I „ji---ifc-i3ifc+i...i„ A . . „. 

tf ^^3l---3k-l^3k+l---3n ^ ^^3l---3k-l33k+l---3n 

has negative eigenvalues. 
Also Corollary Estill holds. 

Corollary 2 A sufficient condition for PPT is that all terms of indexes '2' of the tensor g>^---^" are 0. 

Similarly, if 7r£-) is a permutation of 1, ... ,n, the necessary condition of Proposition |3| generalizes to the sequence of 
inequalities |27| 



tr 



^ tr (pi^^^j^^^^J ^ . . . ^ tr (p2) V permutations tt{-) (22) 
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III. EXAMPLES 

A. Partial quadratic Casimir invariants and classical correlations: an example 

The aim of this example is to show that given a 2-qubit density p, the condition tr (p^) = tr (p^) tr (pg) is not 
a sufficient condition for uncorrelation. Consider the bipartite density: p = p™Aoo + q'^^Aqi + q'^^Aq^ + p^°Aio + 
Q^^Aso + g^^Ai3 + g'^^Asi. Once p"^, p"'^, and are fixed, then g\, g\, g^g and g^g are determined. Choose g^^ 
and g^^ such that 



QaQ% ^ 
QaQb ^ 



but such that 



Then tr {p\) tr (p^) = tr (p^) but C 7^ implying that this state is correlated. Since g^^ = g^"^ 
the correlation is classic. 



0, from Corollary n 



B. Bell state 



The Bell state 1-0) — (|0) + |f)) has density operator given by 



PBell = 2 



1 r 
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In the tensor of coherences, it corresponds to 



PBcu = g' 'Ajk = -Aoo 



Since g- 



y{jk} 7^ {00}, this state has both reduced densities pA and pB that are completely mixed 



tr (p^) — tr (p^) = i. Hence the Bloch vectors gj^ and gs live on "spheres" of radius and in 
J, fc = 1, . . . , 3. tr (pIch) = 1 implies that the entangled state instead is pure. 



Werner states 



Consider the one-parameter family of states 



Pwcr(a;) 



" 1-x 

4 








" 





1+x 
4 


X 

2 








X 

2 


1+x 
4 














1-x 
4 . 



(23) 



which is known to be maximally entangled when x = \ and uncorrelated when a; = 0. Using the computations of 
Appendix 1X1 we get the following nonnuU components: g'^'^ = i, g^^ = — |-, g^^ = — f and g'^^ — 



'|. Hence 



1 X 
PWer{x) = - Aoo ~ 2^^^ 



-A 



22 



-A 



33 



(24) 



For a; = 0, tr (pwgi(0)) — j and, for x = 1, tr (pwer(l)) = 1- From (Pljl . the reduced density operators pA and ps are 
both completely mixed states Va; S [0, 1], since g^'^ = g'^'^ = y{jk} ^ {00}, so, once again, the homogeneous part 
of pwor contains only "pure" correlations: C — pwcr — ^Aqq- 

For this highly symmetric density, the explicit computation of a linear combination of tensor products which is 
separable in the whole domain ^ x ^ i is rather easy with the insight given by the tensor of coherences. For 
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example, we can construct 2-qubit terms g^^ such that the corresponding reduced terms g^'^ and g'^^ are identically 
in the following way. Choose the 12 one-qubit densities: 



PA,1 £'"Ao 


+ gA,i^i 


PA, 2 = fil^Ao + ^^^2^2 


PA,3 




+ ^?A,3A3 






PA, 5 = e^Ao - 6'i^2A2 


PA, 6 = 


g°Xo 


^ eA,3A3 


pb,i ^ g°>^o 


+ eki^i 


PB,2 = ^"Ao + 6l|_2A2 


PB,3 = 


g'Xo 


+ f'l,3A3 


Pea = g^Xo 




PB.h = f^Ao - 6l|_2A2 


PBS = 


g'Xo 


- 6'S,3A3 


equal weights = ... 


vf' = i the 

6 ' 


resulting density p — w^pA p 


<^ Pb,p 


is 




1 , 

P = 2^00 + 


1 1 1 A 


11 + 2£'A,2£'S,2A22 + g Pa, 3^3, 3^33 







(25) 






11 2 2 

Sa.I gfl.l ~gA,2 6B,2 

6 





1 1 I 2 2 
Sa.iBb,1+Ba,2Bb,2 



Sa,iSb,1+8a,2Sb,2 

6 





Ba,3Bb,3 

6 





1 1 _ 2 2 
gA,lgB,l gA.2gB.2 

6 




3 3 
1 I Sa, 35-8.3 

4 6 



While the "multiplicative" part of p is nontrivial, the corresponding reduced densities are always completely mixed 
because of the cancellations occurring. The eigenvalues of H25|) are 



1 1 2 2 3 3 \ 

"£'a,i£'_b,i ~ Qa,2Qb,2 ~ 6a,36b,3) 



' 4 



1 1 2 2 3 3 \ 

"fA.lCs.l + Pa, 2^5, 2 + £'A,3£'B,3j 



\ + \ ((?A,ie'ki - 6'A,2e's,2 + £'A,3£'I,3) ' ^ + ^ {qa,iQba + gA^gB^ - (?A,3£'I,3) > I 



(26) 



All four eigenvalues are ^ for all admissible gA.p, gB,p ^ In particular p of (|25|) is equal to pwcr provided one 

g's 3 = — Qa 3- In tl^^^ case, the eigenvalues never 



— ^3 



and 



-Pa, 



,2 ^ ^'A,2 
2 / 



chooses g\ -^ = 2 = e'A,3 ""-^ t^i3,i 

reach the critical points for a density (0 and +1). If x = — |(gi\ i) — ^{g\ 1)^ then the allowed parameter span 
is a; G [0, as is known for a separable Werner state. In fact, a; ^ by definition (since x — ^{g\ ^ 0), and 
X = I (i) when -j^ reaches i"^- It is interesting to see what happens to the single terms of the convex sum. For 
example 

PA,i«'Ps,i = (e^Ao + eJi^Ai) (g) (e°Ao + £»kiAi) 

= g°g° AoQ + g^gs.i^oi + Qa iQ^^io + Qa ifki^ii 



1 , 
- 2^00 



^ki 



Aoi 



e'A.i 



Aio + £'A,i£'kiAii 



^/2 ' V2 

is a well-defined density within the same boundary as pa,i, Pb,i- Instead if one considers the sum of two terms which 
cancels the one-qubit terms (for example "1" and "4" with equal weights — w^) 

2w^ (eV"Aoo + fi'kieki^ii) 
1 . 3 



W'Va,1 PS,1 + W'^PAA PBA 



= 2w^ i -Aoo 



^xAii ) = 2w^ (pab) 



then pab is a well-defined density when j 



^ 1 i.e., in the larger interval ^ x ^ The same bound 



IS 



found on the other pairwise sums. This is also the bound found by the trace-square inequalities (|15|l . 



Notice that choosing g\ ^ 



s\,2 



g\ 3 but instead g^g ^ 



-Qa2 and g%^ = g\^ and x 



{g\ iPs 1) = ^%{g\ 26's 2)' in the same range of x one gets the different family: 



PWc 



l+x 
4 








4 









4 





2 




l+x 
4 . 



which can also be obtained from pwcr by local unitary transformations and which leads to pBcii considered before for 
X = 1. 
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D. A tripartite family 



-,000 a 1 „ 
5 Aooo + Q 


111 t 


111 — 


1 A 


000 + 


xAiii 


or written 




- 1 

2 V2 





n 





n 


n 





a; 







1 

2 V2 


n 

u 



u 


n 


n 


X 













1 

2 \/2 
















1 











1 

2 \/2 


X 






















X 


1 

2 \/I 



















X 








1 

2 V2 













X 














1 

2 V2 







X 




















1 

2 V2- 



(27) 



^ a; ^ 



has tr (p^) = i + a;^ and eigenvalues || i ^^a;| each with multipUcity 4. p is a density for — ^ ^ 2^2' 

From Corollaryl^l this state is certainly PPT. Furthermore, the necessary conditions (|22|l are satisfied as well as the 
realignment criterion |^ 0, IT^ . Still one may wonder if p is separable and if an explicit mixture can be found and 
how. For X = (and near it) it is obviously separable. Once we have an explicit linear combination for p, then 
this will provide an estimate of its separability interval, valid at least with respect to the mixture used. All 1-qubit 
and 2-qubit reduced densities are in the completely random state. The trick of pairwise canceling reduced densities 
does not generalize in a straightforward manner to multipartite qubits, but some variants of it can still be used. For 
example the following 9 triples of operators 



PAA = g"^o PBS = fi^Ao + QbiXi 

PA,2 = e^Ao pb,2 = ^"Ao - Qb lAl 



PCS = fil^Ao + Qc^iXl 
PC,2 



— f?°Ao — ^?c',lAl 



PA,3 = Q°Xo + QA^a^l 
PA,4 = e^Ao - gA,3Xl 



PS, 3 = e^Ao Pas = £'°Ao + PasAi 

PB,4 = e^Ao pc,4 = e^Ao - Qc^s^l 



PA^b = 6'"Ao + e^gAi 

PAfi fii^Ao - PakAi 



Pb,5 = e'"Ao + gs^gAi 
PB.6 = ei^Ao - 5A1 



Pes, = ^"Ao 



PA,7 = e^Ao + PajAi pb.7 = e^Ao + gsjXi pc.7 = Q^Xq + Qc^rXi 
PAS = g°Xo - Qa.tXi Pbs = ^"Ao - QbjXi pc^ = Q^Xq 

Pa,9 = e°Xo pb,9 = e^Xo pao = e°Xo - Qc^jXi 

with the choices of weights = , = w"^ , = and w'' — ufi = uP such that + . . . + = 1 gives 

^ " 271^°"° ^ 71 '^^^^^BsQcA + Qb, 700,7) Aon + (2w^£'A,3ec,3 + ^«^£'a,7£'c,7) Aioi 
+ \/2 (w^fii^^se'ks + ■f^'^fi'^rek?) Alio + Qa,7Qb, 760, 7^111 
In order to cancel the terms along Aqu, Aiqi, Aho, we fix the parameters of the reduced operators as follows 



£•3.1 



' (x) 3 [w'^) 3 



2u)i 



/ (x) 3 (w'^) 3 



2^3 ' 



1 _J(x)^uF)^ 1 

£'A,5 — V ~ ; 6'S,5 



"Pa, 5' £'yl,7 ^ £'S,7 — BC,7 — \ 
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hence reobtaining (|27|l for all x. What is left to do is to control for which x there exist admissible weights such 
that all the reduced operators are well-defined densities, i.e., Iq]^^] ^ K = A,B,C, r = 1,3,5,7. Squaring the 
relations above, we obtain that x has to fulfill the 4 constraints 

-jr2^-^^2 (28) 

{{xfw'')^ ^w^ , {{x)^w'y {{xfyj'^y (29) 



One can rewrite the last three inequalities as 



, {W^f («;3)3 (u;5)3- 



Using the rule of thumb of selecting weights so that all three reduced densities constraints H3U|I become active simulta- 
neously as we vary x, = w'^ = one gets w"^ = As we vary in ^ ^ | — e , e > 0, the inequalities 
(|28|l and (|30|l give a trade off. A quick numerical search shows that —0.050 ^ x ^ 0.050 is the largest common range 
satisfying the constraints one can achieve, corresponding to = 0.0715. In fact, in correspondence of these values 
we have that all 1- and 2- qubit densities have trace squares between 0.99 and 1 (while tr (p^) = 0.125). 

One may wonder if the bound found depends on the linear combination chosen or less. In order to verify this, it is 
possible to use the alternative family of densities: 

PA,\ fii^Ao pb.\ = ^"Ao PCS = q">^q - Qcs^i 

PAa = £'°Ao + 0X2^1 Pb,2 = ^"Ao -I- Ps^aAi Pea = £'"Ao 
PA,?, = £'°Ao - 0X2^1 PB,3 = ^"Ao - 6'b,2Ai Pas = £'"Ao 

PAA = e°Ao + Qaa^I PBA = e^^O + 6'i3,4Al PCA = e°>^0 + QCA^I 

PA,5 = e°Ao + g\Ai pb,5 = ^"Ao + Qb 4A1 Pc,5 = g°>^o - QCA^I 



corresponding to 



+ 2w**£'A,46'S,4£'C,4Alll 

Choosing = w"^, = w'^ and 



1 1 4 11 /~ /91 1 4-1 1\ 

^000 + 7. + W Qc,i Aool + V2 (w Qa,2Sb,2 + ^ SA,iQB,i) ^11 



Be, I 



QA,2 



Pa, 4 — &B,A — Oca 




this convex combination has components that are all densities when the three inequalities are satisfied: 



%y/2{w^f 8\/2(w4)2' V 2u;4 ^ " V 2^4 ' 

A parametric search gives that x is separable in —0.1166 ^ a; ^ 0.1166 in correspondence for example of = 0.33, 
vY = 0.17 and = 0.165. Hence just like the mixture representing a density is nonunique, also the separability 
range changes with the convex combination chosen. No conclusion is drawn about separability in 0.1166 ^ ^ 
although this density is probably separable on the whole range. 
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E. Another tripartite biseparable family 

Consider the 3-qubit state p — q''^'- Kjki where g^*^^ = and 



Expanding exphcitly 



P = 



j031 _ 










= £)303 






f?330^ 


^>33l = 


jOll _ 










= £)301 


= —X 








pkl ^ 


otherwise. 


















r\/2+12 

4 


X 


a; 


— X 


X 


— X 


— X 


X 






X 


4 


— X 


X 


— X 


—3 X 


X 


— X 






a; 


—a; 


\/2-4 2; 
4 


—3 X 


— X 


X 


X 


— X 


1 




—a: 


a; 


—3 X 


V2-4£C 
4 


X 


— X 


— X 


X 






X 


— X 


— X 


X 


\/2-4 3; 
4 


X 


—3 X 


— X 






— X 


—3 X 


X 


— X 


X 


\/2-4 2; 
4 


— X 


X 






— X 


a; 


X 


— X 


—3 X 


— X 


\/2-4 2; 
4 


X 






X 


—a; 


— X 


X 


— X 


X 


X 


\/2+12a: 
4 -1 



(31) 



The trace norm for this density is tr (p^) = | + 16x^ and the eigenvalues {| ± V2a;} each with muhiphcity 4. Hence p 
is a well defined density matrix for ^ 2; ^ and it is mixed in the entire interval (completely mixed for x = 0). 

When X ~ we recover the example of Eq. (6) of 4] puPB = \ (^Is — j with {ipj) an Unextendible 

Product Basis (UPB) for 3-qubit states: = |01+), |1 + 0), | + 01), | ) (where |±) = (|0) ± |1))). puPB is 

known to be entangled and PPT. Again, this last fact is simply verified by noticing that no index "2" appear in l|31(l . 
hence Corollary |21 holds for all x. 

Notice that for this example the necessary conditions H22I) are fulfilled for all x, but that p is separable only in some 
proper subinterval [— Xc, Xc] around the complete mixing state. 

Also in this case it is possible to obtain a linear convex combination of tensor products of trace 1 Hermitian matrices 
corresponding to p. This allows to estimate the separability interval with respect to the particular mixture chosen by 
imposing that each matrix of each tensor product be a well-defined density. One notices first that 1-qubit reduced 
densities are all completely mixed (g"*" — g'^'^'^ = = 0). Hence for each nonzero 2-qubit term a cancellation 
on the corresponding reduced terms must occur, as in the Werner state parameterization. On the contrary, 2-qubit 
reduced densities are not completely random and it is straightforward to observe that all of the nonzero 2-qubit terms 
correspond to "reductions" of the 4 nonzero 3-qubit terms: for example g^^^ 7^ implies g''^^^ ^ 7^ g^^'^ 7^ 0. 
If we use the same scheme of Section IIII Dl to produce the 3 party terms, then we obtain a convex combination of 
36 triples of two level densities which, with the naive choice of all equal weights, guarantees separability only in a 
tiny interval: —0.006 ^ x ^ 0.006 (the realignment criterion detects entanglement only close to -^^)- The explicit 
expression of the components is available upon request. Notice that the linear combination obtained connects puPB 
with g^^^Aooo- From the classification of [l|, it is known that the biseparable entangled regions have a border in 
common with the set of separable states. Since this last set is convex, the linear combination has a unique crossing 
point Xc' (plus, specularly, —Xc')- 

Managing densities and weights becomes rapidly cumbersome, but one can expect the problem to admit an algo- 
rithmic formulation via (tractable) convex optimization schemes. 

F. Building NPT entangled states 

NPT entanglement (i.e., entanglement detectable by means of the PPT test) is obviously the simplest to manipulate 
for the purposes of building entanglement between subsystems. There are by now standard methods to do that, 
typically starting from particular ket states or taking linear combinations of known singlets with product states. The 
aim of this example is just to show that such a construction is also easy in the formalism we are proposing, and the 
states one obtains are not the standard one normally considered in the literature (GHZ, W states, etc.) 30]. For 
example, assume that one wants to have three qubits with A and B NPT entangled, C satisfying the PPT property. 
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J7 
212 



(a) (b) (c) 

FIG. 1: {x,y) parameter space for the example of Section llll Fl In (a) and (b) the triangle represents the admissible {x,y) 
pairs for pab-c and P^b-c respectively. Hence {x,y) belonging to the complement of the rhomb in the triangle of (c) gives 
NPT entangled densities. 



One possible solution is the following state 

Pab-c = + xKi22 + XA212 + 2/A330 (32) 

For pab-c we have tr {p\b_c) — |+2a;^+j/^ and eigenvalues | (l — 2\/2j/) of multiplicity 4 and | (l + 2\/2 {y ± 2a;)) 
each of multiphcity 2. p^^ and p^^ instead have eigenvalues | (l + 2\/2j/) with multiplicity 4 and | (l — 2%/^ {y ± 2a;)) 

each of multiplicity 2. pab-c is a state in the triangle of FigureHla); while p^b-c (*-"^ Pab-c) positive eigenvalues 
in the flipped triangle of Figure Q] (b). Hence the area outside the rhomb of Figure Q (c) gives the pairs (x, y) leading 
to a NPT entangled state. Notice that dropping any of the three terms in (|32|) the density becomes PPT. 

IV. CONCLUSION 

For multiparty qubit densities, the tensorial formalism proposed, although intuitively simple and certainly not so 
original, has several advantages which are now summarized. 

• First and foremost: it provides a lot of inside in the correlation pattern between subsystems and between 
"groups" of subsystems, as it is known from the literature ,l^,2m. Hence it may be useful in the construction 
of quantum networks. 

• ft allows to explain the concept of reduced density in terms of the afhne parameterization and the tensor itself 
is given by all "degrees of reduction" down to the single qubit densities in an unambiguous way. 

• It also allows to give a geometric picture of the density in terms of juxtaposition of affinc Bloch spheres. Thanks 
to the afhne parameterization, also reduced densities enter into the picture. For example the degree of purity 
of the n-party densities is given in terms of the radius of the sphere of the corresponding dimension via the 
standard trace norm. Such norm itself can be intended as a Tsallis entropy for a particular choice of the index 
or as a quadratic Casimir invariant for the set of density operators. 

• It allows to construct a density with the desired correlations between its subsystems. In particular also quantum 
correlations of the NPT type. 

• It provides an intuitive interpretation of both NPT and bound entanglement: they correspond to linear combi- 
nations of tensor products in which at least one of the factors in one of the products is not a well-defined density 
in the sense that its Bloch vector is too big in norm. 
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In conclusion, we hope that the formahsm put forward in this paper may help not only in the understanding of 
quantum correlations in multipartite systems, but also in their concrete and systematic engineering. 



APPENDIX A: TENSOR OF COHERENCES FOR TWO QUBIT STATES 

If {p)pq p, g = 1, . . .4, are the elements of the density operator p of Section FlI Bl the tensor of coherences = 

tr (pAjfe), j, fc = 0, . . . , 3, corresponding to it has components: 

= ^ ((P)ll + (P)22 + (P)33 + (P)44) 

= Re[(p)l2] + Re[(p)34] 

= -Im[(p)i2]-Im[(p)34] 

= ^((p)ll-(/5)22 + (p)33-(p)44) 

= Re[(p)i3] + Re[(p)24] 

= -Im[(p)i3] - Im[(p)24] 

= 2 + (^^22 - (p)33 - (P)44) 

= Re[(p)i4] + Re[(p)23] 

= -Im[(p)i4]+Im[(p)23] 

= Rc[(p)i3]-Re[(p)24] 

= -Im[(p)i4]-Im[(p)23] 

= -Rc[(p)i4]+Re[(p)23] 

p23 ^ -Im[(p)i3] +Im[(p)24] 

= Re[(p)i2] - Re[(p)34] 

= -Im[(p)i2] +Im[(p)34] 

Q^^ = ^ ((P)ll - (P)22 - (P)33 + (P)44) 
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